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This paper evolved from a long series of discussions between the two 
authors, going back to around 1980, on the problems of making effective 
computations in algebraic geometry, and it took more definite shape in a 
survey talk given by the second author at a conference on Computer Algebra 
in 1984. The goal at that time was to bring together the perspectives of 
theoretical computer scientists and of working algebraic geometers, while 
laying out what we considered to be the main computational problems and 
bounds on their complexity. Only part of the talk was written down and since 
that time there has been a good deal of progress. However, the material that 
was written up may still serve as a useful introduction to some of the ideas 
and estimates used in this field (at least the editors of this volume think so), 
even though most of the results included here are either published elsewhere, 
or exist as "folk-theorems" by now. 

The article has four sections. The first two parts are concerned with 
the theory of Grobner bases; their construction provides the foundation for 
most computations, and their complexity dominates the complexity of most 
techniques in this area. The first part introduces Grobner bases from a 
geometric point of view, relating them to a number of ideas which we take 
up in more detail in subsequent sections. The second part develops the 
theory of Grobner bases more carefully, from an algebraic point of view. 
It could be read independently, and requires less background. The third 
part is an investigation into bounds in algebraic geometry of relevance to 
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these computations. We focus on the regularity of an algebraic variety (see 



Definition p.2|) , which, beyond its intrinsic interest to algebraic geometers, 
has emerged as a measure of the complexity of computing Grobner bases 
pS87a|| , |pS87b|| , PS88|| ). A principal result in this part is a bound on 



see 



the regularity of any smooth variety by the second author: Theorem |3.12| (b). 



This bound has stimulated subsequent work, and has now been generalized 
by IP3EL91 ]. Another result of interest is Proposition |3.13 , which elucidates 



the scheme structure of the ideal membership problem. The fourth part 
is a short discussion of work on algorithms for performing some other key 
operations on varieties, some open problems about these operations and some 
general ideas about what works and what doesn't, reflecting the prejudices 
of the authors. 

One of the difficulties in surveying this area of research is that mathemati- 
cians from so many specialties have gotten involved, and they tend both to 
publish in their own specialized journals and to have specific agendas corre- 
sponding to their area. Thus one group of researchers, the working algebraic 
geometers, are much more interested in actually computing examples than 
in worst-case complexity bounds. This group, including the first author, has 
put a great deal of work into building a functioning system, Macaulay, based 
on Grobner bases, which has solved many problems and provided many ex- 
amples to the algebraic geometry community |[BS92a|]. Another group comes 



from theoretical computer science and is much more interested in theoretical 
bounds than practical systems (cf. the provocative comments in Lenstra's 
survey |[Len92|| ). It seems to us that more communication would be very 



helpful: On the one hand, the working algebraic geometer knows lots of facts 
about varieties that can be very relevant to finding fast algorithms. Con- 
versely asymptotic and/or worst-case performance bounds are sometimes, 
at least, important indicators of real-time performance. These theoretical 
bounds may also reveal important distinctions between classes of procedures, 
and may pose new and deep problems in algebraic geometry. Thus we will 
see in Section |^ how regularity estimates fiesh out a picture explaining why 
Grobner basis computations can have such explosive worst case behavior, yet 
be so useful for the kinds of problems typically posed by mathematicians. Fi- 
nally, to make this article more useful in bridging this gap, we have tried to 
include a substantial number of references in our discussions below. 
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1 A Geometric Introduction 



Let X be a subvariety or a subscheme of projective n-space P", over a field 
k. Let T he a, vector bundle or a coherent sheaf supported on X. We would 
like to be able to manipulate such objects by computer. From algebra we get 
finite descriptions, amenable to such manipulations: Let S = k[xo, . . . ,Xn] 
be the homogeneous coordinate ring of P". Then X can be taken to be the 
subscheme defined by a homogenous ideal I (Z S, and T can be taken to be 
the sheaf associated to a finitely generated ^'-module M. We can represent / 
by a list of generators {fi, . . . , f^), and M by a presentation matrix F, where 

Ml Mo — ^ M — ^ 

presents M as a quotient of finitely generated free 5*- modules Mq, Mi. We 
concentrate on the case of an ideal /; by working with the submodule J — 
lm{F) C Mo, the module case follows similarly. 

The heart of most computations in this setting is a deformation of the in- 
put data to simpler data, combinatorial in nature: We want to move through 
a family of linear transformations of P'* so that in the limit our objects are 
described by monomials. Via this family, wc hope to pull back as much 
information as possible to the original objects of study. 

Choose a one-parameter subgroup X{t) C GL{n + 1) of the diagonal form 

where W = {wq, . . . , w„) is a vector of integer weights. For each t 7^ 0, X{t) 
acts on X via a linear change of coordinates of P", to yield the subscheme 
Xt = X{t)X ^ X. The hmit 

Xo = \imXt 

is usually a simpler object, preferable to X for many computational purposes. 

Even if we start out by restricting X to be a subvariety rather than a sub- 
scheme of P", it does not suffice to take the limit Xq set-theoretically; often 



\{t) = 
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all we will get pointwise in the limit is a linear subspace L C P", reflecting 
little besides the dimension of the original variety X. By instead allowing 
this limit to acquire embedded components and a nonreduced structure, we 
can obtain an Xq which reflects much more closely the character of X itself. 

We compute explicitly with the generators fi, . . . , of /: Let A act on 
S by mapping each Xi to t^^xf, A maps each monomial x"^ = Xg" ■ ■ ■ a;^" 
to t^-^x^ = t'"o^o+-+vjnanx'^o ...x^". If / = ax^ + 6x^ + . . ., then A/ = 
at'^-^x^ + 6t^-^x^ + . . .. We take the projective limit in(/) = limt^o A/ 
by collecting the terms of A/ involving the least power oit; in(/) is then the 
sum of the terms ax^ of / so • A is minimal. For a given / and most 
choices of A, in(/) consists of a single term. 

The limit Xq we want is deflned with all its scheme structure by the 
ideal in(/) = lim^^^o A/, generated by the set { in(/) | / G / }. For a given 
/ and most choices of A, in(/) is generated by monomials. Unfortunately, 
this deflnition is computationally unworkable because / is an inflnite set, and 
in(/) need not equal (in(/i), . . . , in(/r)) for a given set of generators /i, . . . , 
of /. To understand how to compute in(7), we need to look more closely at 
the family of schemes Xt deflned by A. 

Let S[t] be the polynomial ring k[xo, . . . , Xn, t]; we view S[t] as the coor- 
dinate ring of a one-parameter family of projective spaces P" over the affine 
line with parameter t. For each generator fj of /, rescale Xfj so the lowest 
power of t has exponent zero: Let Qj — f^Xfj, where £ — W ■ A is the least 
exponent of t in Xfj. Then fj = gj\t=i and in(/j) = gj\t=o- Now, let J C S[t] 
be the ideal generated by {gi, . . . , gr); J deflnes a family Y over whose 
central flber is cut out by (in(/i), . . . , in(/r)). 

What is wrong with the family Y7 Y can have extra components over 
t — 0, which bear no relation to its hmiting behavior as i — > 0. Just as 
the set-theoretic limit limt_^o can be too small (we need the nonreduced 
structure), this algebraically deflned limit can be too big; the natural limit 
lies somewhere in between. 

The notion of a flat family captures exactly what we are looking for here. 
For example, if Y is flat, then there are no extra components over t = 0. 
While the various technical deflnitions of flatness can look daunting to the 
newcomer, intuitively flatness captures exactly the idea that every flber of a 
family is the natural scheme-theoretic continuation of its neighboring flbers. 
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In our setting, all the Xt are isomorphic for t 7^ 0, so we only need 
to consider flatness in a neighborhood of t = 0. Artin ||Art76 | gives a cri- 



terion for flatness applicable here: The syzygies of gi, . . . ,gr are the rela- 
tions higi + . . . + hrgr = for hi, . . . ,hr G S[t]. Syzygies correspond to 
elements {hi,...,hr) of the S'[t]-module S[tY; the set of all syzygies is a 
submodule of 5'[t]''. F is a fiat family at t = if and only if the restric- 
tions {hi\t=o, . . . , hr\t=o) of these syzygies to the central fiber generate the 
S'-module of syzygies of gi\t=o, • • • , 5'r|t=o- 

When gi \t=o, ■ ■ ■ ,gr \t=o are single terms, their syzygies take on a very 
simple form: The module of syzygies of two terms ax^, is generated 
by the syzygy bx'"{ax^) — ax^{bx^) = 0, where = x'^x'^ = x^x^ is 
the least common multiple of x"^ and x^. The module of syzygies of r such 
terms is generated (usually not minimally) by the syzygies on all such pairs. 

We want to lift these syzygies to syzygies of gi, . . . , gr, working modulo 
increasing powers of t until each syzygy lifts completely. Whenever we get 
stuck, we will find ourselves staring at a new polynomial gr+i so t^gr+i G J 
for some i > 0. Including g^+i in the definition of a new J' D J has no 
effect on the family defined away from t = 0, but will cut away unwanted 
portions of the central fiber; what we are doing is removing t-torsion. By 
iterating this process until every syzygy lifts, we obtain explicit generators 
gi, . . . ,gr, gr+i, gs for a. fiat family describing the degeneration of X = Xi 
to a good central fiber Xq. The corresponding generators gi\t=i, • • • , gs\t=i of 
/ are known as a Grobner basis for /. 

This process is best illustrated by an example. Let S = k[w,x,y, z] be 
the coordinate ring of P'^, and let I = (/i, /2, /s) C S for 

h = w'^ - xy, f2 = wy- xz, f^ = wz- y'^. 

I defines a twisted cubic curve X C P^; X is the image of the map (r, s) 1— >■ 
(r^s, r^, rs^, s^). Let 



t 



-16 



If w°'x^y^z'^ is a monomial of degree < 4, then A • w'^x^y'^z'^ = t ^w°'x^y'^z'^ 
where ^ = 16a + 46 + c. Thus, sorting the monomials of 5* of each degree 
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< 4 by increasing powers of t with respect to the action of A is equivalent to 
sorting the monomials of each degree in lexicographic order. 

We have 

93 = ^ wz- t'V. 

The module of syzygies on w"^, wy, wz is generated by the three possible 
pairwise syzygies; we start with the syzygy y{w^) —w{wy) — 0. Substituting 
Qi, g2 for the lead terms w'^, wy we get 

y{w^ — t^'^xy) — w(wy — t^^xz) = t^^wxz — t^'^xy^ 

which is a multiple t^^x of g^. Thus, the syzygy 

ygi - W92 - t^^xg^ = 

of fl'ij ^72! 93 restricts to the monomial syzygy y(w^) — w{wy) — when we 
substitute t = 0, as desired. 

Similarly, the syzygy 

zgi - t^Sg2 -wgs^O 

restricts to the monomial syzygy z{w'^) — w{wz) = 0. When we attempt to 
lift z{wy) — y{wz) = 0, however, we find that 

z{wy - t^^xz) - y{wz - t^^y"^) = -t^^xz^ + t^^y^. 

xz^ is not a multiple of w^, wy^ or wz^ so we cannot continue; J — {gi, g2, gs) 
does not define a fiat family. Setting t — 1, the troublesome remainder is 
—xz'^ + y^. Making this monic, let — xz^ — y^; f4 e I and 

g^ = t^A/4 = xz'^ - ty^. 

Adjoin g4 to the ideal J, redefining the family Y. Now, 

zg2 - y93 + 1^^94 = 
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restricts to z(wy) — y{wz) = as desired. 

The module of syzygies of w"^, wy, wz, and xz^ is generated by the 
pairwise syzygies we have already considered, and by the syzygy xz{wz) — 
w{xz^) — 0, which is the restriction of 

-ty'^g2 + xzgs - wg^ = 0. 

Thus, J — {gi, g2, gs, g4:) defines a flat family Y, and 

— xy, wy — xz, wz — y^, x^ — y^ 

is a Grobner basis for /. The limit Xq is cut out by the monomial ideal 
in(/) = {w'^,wy,wz,xz'^), which we shall see shares many properties with 
the original ideal /. Note that xz"^ — y^ = defines the projection of X to 
the plane in x, y, and z. 

The scheme structure of Xq is closely related to the combinatorial struc- 
ture of the monomial /c-basis for S'/in(/): For each degree d in our example, 
the monomials not belonging to in(/) consist of three sets {x'^, x'^^^y, . . . , y*^}, 
{x'^~^z, x'^~'^yz, . . . , y'^~^z}, {y'^, y'^~^z, . . . , z'^}, and a lone extra monomial 
x'^~^w. The first two sets correspond to a double line supported on -u; = 2; = 0, 
the third set to the line w = x = 0, and the extra monomial to an embedded 
point supported at w — y — z — 0. Together, this describes the scheme 
structure of Xq. The first two sets consist of 1 and d monomials, respec- 
tively; the third set adds d — 1 new monomials, and overlaps two monomials 
we have already seen. With the extra monomial, we count 3d-\-l monomials 
in each degree, which agrees with the dimensions of the graded pieces of S/ 1. 
The embedded point is crucial; it makes this count come out right, and it 
alone keeps Xq nonplanar like X. 

The new monomial generator xz^ of in(7) excludes the line w — y — 
from Xq] combinatorially, it excludes all but three monomials of the set 
{x'^ , x'^~^ z , . . . , z'^} from the monomial /c-basis for each degree of the quo- 
tient 5'/in(/). We can see that this line is unwanted as follows: Away from 
t — 0, Y is parametrized by (r, 1— > (t^^r^s, t'^r^, trs^, s^, t). Thus, fix- 
ing r and s, the curve {r,ts,t) 1— > (i^V^s, iV^, t^rs^, t^s^, t), with projec- 
tive limit (0,0,r, s,0) as t — >■ 0. Similarly, the curve {r,t^s,t^) has as its 
limit (0, r^, s^, 0, 0). These calculations show that the lines w = z = Q and 
w = X = Q indeed belong set-theoretically to the limit Xq. We can find no 
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such curve whose hmit is a general point on the hne w = y = 0, for (r, t^s, t^) 
doesn't work. Thus, the hne w = y = sticks out of the good total space Y. 

One usually computes Grobner bases by working directly in the ring S, 
dispensing with the parameter t. The one-parameter subgroup A is replaced 
by a total order on the monomials of each degree, satisfying the multiplicative 
property x"^ > x-^ =^ x'-^x"^ > x'-^x'^ for all x'-^. In fact, for our purposes 
these are equivalent concepts: The weight vector W associated with A induces 
the order x"^ > x"^ <^==^ W-A < W-B, which is a total multiplicative order in 
low degrees as long as no two monomials have the same weight. Conversely, 
given any multiplicative order and a degree bound d, one can find many A 



which induce this order on all monomials of degree < d. See Pay82|| , |[Rob85 
for characterizations of such orders. 

We shall be particularly interested in two multiplicative orders, the lex- 
icographic order used in our example, and the reverse lexicographic order. 
The lexicographic order simply expands out the monomials of each degree 
into words, and sorts them alphabetically, i.e. x"^ > x"^ iff the first nonzero 
entry in A — S is positive. The reverse lexicographic order pushes highest 
powers of x„ in any expression back to the end, then within these groups 
pushes highest powers of Xn-i to the end, etc., i.e. x"^ > x^ iff the last 
nonzero entry of A — i? is negative. 

What do these orders mean geometrically? The dominant effect of the 
lexicographic order is a projection from P" to P"~^, eliminating xq. A second 
order effect is a projection to P"~^, and so forth. We could compute the 
deformation from X to Xq with respect to the lexicographic order in stages 
carrying out these projections, first applying a A with W = (—1,0, . . . ,0), 
then with W = {—1, —1,0, . . . ,0), etc. Alternatively, for monomials of each 
degree < d, we can apply the single A with W = (— . . . , — rf, — 1, 0), 
generalizing the A used in our example. Use of the lexicographic order tends 
to muck up the family Y more than necessary in most applications, because 
projections tend to complicate varieties. 

For the reverse lexicographic order, the dominant effect is a projection 
of P" down to the last coordinate point (0, . . . , 0, 1). As a secondary effect, 
this order projects down to the last coordinate line, and so forth. In other 
words, this order first tries to make X into a cone over the last coordinate 
point, and only then tries to squash the result down to or cone it over the 
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last coordinate line, etc. For monomials of each degree < d, this can be 
realized by applying A with W = (0, 1, c?, . . . , rf""^). Like such cones, the 
reverse lexicographic order enjoys special properties with respect to taking 
linear sections of X or Xq by intersection with the spaces defined by the last 



variable(s) (see |pS87a|| ). The preferred status of the reverse lexicographic 
order can be attributed to this relationship, because generic linear sections 
do not complicate varieties. 

For example, if we take X to be three general points in P^, then using 
the lexicographic order Xq becomes a triple point on a line, because the first 
order effect is the projection of the three points to a line, and the second 
order limiting process keeps the points within this line. By contrast, if we 
use the reverse lexicographic order then Xq becomes the complete first order 
neighborhood of a point (a point doubled in all directions). This is because 
the first order limiting process brings the three points together from distinct 
directions, tracing out a cone over the three points. The first order neigh- 
borhood of the vertex in this cone has multiplicity 3, and is the same as the 
complete first order neighborhood in the plane of this vertex. 

For those familiar with the theory of valuations in birational geometry 



^S76| , Vol. II, Ch. VI], the lexicographic and reverse lexicographic orders 



have simple interpretations. Recall that if X is a variety of dimension n, and 
F -.X = ZoD ZiD Z2D ...D Zn 

is a flag of subvarieties, codimx(^i) = i, with Zi smooth at the generic point 
of Zj+i, then we can define a rank n valuation vp on X as follows: For each 
i = 1, . . . , — 1, fix /j to be a function on with a 1^^ order zero on Z^. 
Then for any function /, we can define ei = ord^^ (/) , 62 = ord^j ( (///i ^ ) Ui ) , 
etc., and vp{f) = (ei,...,e„) G Z", where the value group Z" is ordered 
lexicographically. The arbitrarily chosen /j are not needed to compare two 
functions /, g: We have vpif) >- vpig) if and only if ordziif / g) > 0, or if this 
order is zero and ordz2((//5') UJ > 0, and so forth. Such a valuation also 
defines an order on each graded piece Sd of the homogeneous coordinate ring: 
take any /o G 5"^ and say f > g if and only if vp{f/fo) y vpig/fo). More 
generally, one may take the Zi to be subvarieties of a variety X' dominating 
X and pull back functions to X' before computing vp- 

The lexicographic order on monomials of each degree of P" is now induced 
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by the flag 

P" D V{xo) D V{xo, Xi)D ...D V{xo, Xn-l). 

For example, the first step in the comparison defining vp^^k^/Jq) >- 
vf{^^ / fo) has the effect of asking if qq — bo > 0. 

The reverse lexicographic order is induced by a flag on a blowup X of 
P": First blow up V{xo, . . . ,Xn-i) and let Ei be the exceptional divisor. 
Next blow up the proper transform of ^(a^o, . . . , Xn-2), and let E2 be this 
exceptional divisor. Iterating, we can define a flag 

X D EiD EinE2D ...D Ein...nEn 

which induces the reverse lexicographic order on monomials in each degree. 
For example, looking at the affine piece of the first blow up obtained by 
substituting xq — x^Xn-i, ■ ■ ■ , = x'n-2^n-i, the power of Xn-i in the 

transform of x"^ is Oq + . . . + cin-i, which is the order of vanishing of this 
monomial on Ei. Thus, the first step in the comparison defining Vi?(x"^//o) >^ 
fir(x^//o) has the effect of asking if ao + . . . + a„_i — bo — ■ ■ ■ — fcn-i > 0, 
which is what we want. 

Taking into account the equivalence between multiplicative orders and 
one-parameter subgroups, the process we have described in S[t] is exactly 
the usual algorithm for computing Grobner bases. It is computationally 

advantageous to set t = 1 and dismiss our extra structure as unnecessary 
scaffolding, but it is conceptually advantageous to treat our viewpoint as 
what is "really" going on; many techniques of algebraic geometry become 
applicable to the family Y, and assist in analyzing the complexity of Grobner 
bases. Moreover, this picture may help guide improvements to the basic 
algorithm. For example, for very large problems, it could be computationally 
more efficient to degenerate to Xo in several stages; this has not been tried 
in practice. 

The coarsest measure of the complexity of a Grobner basis is its maximum 
degree, which is the highest degree of a generator of the ideal in(7) defining 
Xo- This quantity is bounded by the better-behaved regularity of in(/): The 
regularity of an ideal / is the maximum over all i of the degree minus i of any 
minimal i^^ syzygy of /, treating generators as 0*^ syzygies. When / is the 
largest (the saturated) ideal defining a scheme X, we call this the regularity 
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of X. We take up regularity in detail in Section here it suffices to know- 
that regularity is upper semi- continuous on flat families, i.e. the regularity 
can only stay the same or go up at special fibers. 

Let reg(/) denote the regularity of J, and regQ(J) denote the highest 
degree of a generator of /. In our case, t = is the only special fiber, and 
the above says that 

ego(/) < reg(/) < reg(in(/)) > rego(in(/)), 

where rego(/) can be immediately determined from the input data, and 
rego(in(/)) is the degree-complexity of the Grobner basis computation. In 
practice, each of these inequalities are often strict. 

However when k is infinite, then for any set of coordinates for P" chosen 
from a dense open set U C GL{n+l) of possibilities, Galligo ( ||Gal74|| ; see also 



ret 



[ [BS87b| ) has shown that the limiting ideal in(/) takes on a very special form: 
in(/) is invariant under the action of the Borel subgroup of upper triangular 
matrices in GL{n + 1). This imposes strong geometric conditions on Xq. In 
particular, the associated primes of in(J) are also Borel-fixed, so they are all 
of the form (xq, . . . , Xj) for various i. This means that the components of Xq 
are supported on members of a flag. 



In characteristic zero, it is shown in pS87a|| that the regularity of a Borel- 
fixed ideal is exactly the maximum of the degrees of its generators, or in our 
notation, that reg(in(J)) = rego(in(/)) when in(J) is Borel-fixed. Thus, for 
generic coordinates in characteristic zero, the degree-complexity of comput- 
ing Grobner bases breaks down into two effects: the gap regQ(/) < reg(/) 
between the input degrees and the regularity of X, and the gap reg(/) < 
reg(in(/)) allowed by upper-semicontinuity. 

A combination of theoretical results, hunches and experience guides the 
practitioner in assessing the ffist gap; what about the second? Does the 
regularity have to jump at all? One can easily find examples of ideals and 



total orders exhibiting such a jump, but in |BS87a|, it is shown that for the 



reverse lexicographic order, in generic coordinates and any characteristic, 
there is no jump: reg(/) = reg(in(/)), so in characteristic zero we have 

rego(in(J)) = reg(J). 

In this sense, this order is an optimal choice: For the reverse lexicographic 
order, the degree- complexity of a Grobner basis computation is exactly the 
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regularity of the input data. This agrees with experience; computations made 
on the same inputs using the lexicographic order can chmb to much higher 
degrees than the reverse lexicographic order, in practice. 

For many applications, one is free to choose any order, but some prob- 
lems restrict us to using orders satisfying combinatorial properties which the 
reverse lexicographic order fails to satisfy. An example, developed further in 
Section H, is that of eliminating variables, or equivalently, of computing pro- 
jections. To compute the intersection of / with a subring R = k[xi, . . . , Xn], 
it is necessary to use an order which in each degree sorts all monomials not 
in R ahead of any monomial in R. The lexicographic order is an example 
of such an order, for each i simultaneously. This strength comes at a cost; 
we are paying in regularity gaps for properties we may not need in a partic- 
ular problem. An optimal order if you need one specific projection (in the 
same sense as above) is constructed by sorting monomials by total degree 
in the variables to be eliminated, and then breaking ties using the reverse 



lexicographic order. See [pS87b|| for this result, and a generalization to the 
problem of optimally refining any nonstrict order. 

Using this elimination order, one finds that the inherent degree- 
complexity of a computation is given not by the regularity of X itself, but 
rather by the regularity of the flat projection X' of X, which is the central 
fiber of a fiat family which animates the desired projection of X as t ^ 0. 
The jump in regularity between X and X' is unavoidable; by choosing an 
optimal order, we avoid the penalty of a further jump in regularity between 
X' and Xq. 

The regularity of algebraic varieties or schemes X is far from being well 
understood, but there is considerable interest in its study; this computational 
interpretation of regularity as the inherent degree-complexity of an ideal is 
but one more log on the fire. 

From a theoretical computer science perspective, the full complexity of 
computing Grobner bases is determined not merely by the highest degree 
regQ(/) in the basis, but by the total number of arithmetic operations in the 
field k required to compute this basis. This has not been analyzed in general. 



but for 0-dimensional ideals /, Lakshman and hazard ( [ [Lak9 If , ||LL91|| ) have 



shown that the complexity of computing reduced Grobner bases is bounded 
by a polynomial in rf", where d is the maximum degree of the generators. 
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and n is the number of variables. 



2 Grobner Bases 



Let S = k[xo, . . . ,Xn] be a graded polynomial ring over the field k, and let 
/ C S" be a homogeneous ideal. 

Let Sd denote the finite vector space of all homogeneous, degree d poly- 
nomials in S, so S = SqQ) Si Q) ... ® Sd® .. .. Writing I in the same manner 
as / = Jo © /i © . . . © /d © . . . , we have Id C 5*^ for each d. Recall that the 
Hilbert function of / is defined to be the function p{d) = dim{Id), for d > 0. 

A total order > on the monomials of S is said to be multiplicative if 
whenever x"^ > x'^ for two monomials x"^, x"^, then x'^x^ > x'-^x'^ for all 
monomials x*^. This condition insures that if the terms of a polynomial are 
in order with respect to >, then they remain in order after multiplication by 
a monomial. 



Definition 2.1 Let > be a multiplicative order. For a homogeneous polyno- 
mial f = Cix"^! + . . . + CmX"^™ with x"^! > . . . > x"^™, define the initial term 
in(/) to be the lead (that is, the largest) term cix^^ of f . For a homogeneous 
ideal I G S , define the initial ideal in(/) to be the monomial ideal generated 
by the lead terms of all elements of I. 

Note that the definitions of in(/) and in(J) depend on the choice of mul- 
tiplicative order >. See [[BM88|] and [[MR88|| for characterizations of the finite 



set of in(J) realized as the order > varies. 
Fix a multiplicative order > on 5*. 

Proposition 2.2 (Macaulay) / and have the same Hilbert function. 

Proof. ([ [Mac27|| ) The lead terms of Id span m{I)d, because every monomial 
x"^ G in(/) is itself the lead term in(/) of some polynomial f E I: Since 
_ -^c^B some x"^ = in{g) with g E I, we have x"^ = in(/) for 

Choose a fc-basis Bd C Id with distinct lead terms, and let in{Bd) be the 
set of lead terms of Bd; in{Bd) has cardinality p{d) = dim(/rf). Since any 
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element of Id is a linear combination of elements of Bd, any lead term of Id is 
a scalar multiple of an element of in{Bd). Thus, m{Bd) is a basis for in{I)d, 
so p{d) = dim(in(/)(i). | 

One can compute the Hilbert function of / by finding in(/) and applying 
this result; see |PM88|| , §URgT[, and fBS^. 



Corollary 2.3 The monomials of S which don't belong to in(J) form a k- 
basis for S/I. 

Proof. These monomials are linearly independent in S/I, because any linear 
relation among them is a polynomial belonging to /, and all such polynomials 
have lead terms belonging to in(J). These monomials can be seen to span 



S/I by a dimension count, applying Proposition pT2 



Two examples of multiplicative orders are the lexicographic order and the 
reverse lexicographic order, > x^ in the lexicographic order if the first 
nonzero coordinate of A — i? is positive. For example, if S" = k[w,x,y, z], 
then w > X > y > z in Si, and 

w"^ > wx > wy > wz > x^ > xy > xz > y"^ > yz > z^ 

in S\ 

^ reverse lexicographic order if the last nonzero coordinate 

of A — 5 is negative. For example, if S" = k[w, x, y, z], then w > x > y > z 
in Si, and 

w'^ > wx > x^ > wy > xy > y"^ > wz > xz > yz > z"^ 

in S"^. These two orders agree on 5*1, but differ on the monomials of S of 
degree > 1 when n > 2. 

The lexicographic order has the property that for each subring 
k[xi, . . . , Xn] C S and each polynomial / G 5*, / G k[xi, . . . , x„] if and only if 
in(/) G k[xi, . . . , Xn]- The reverse lexicographic order has the property that 
for each / G k[xo, . . . ,Xi], Xi divides / if and only if Xi divides in(/). 

One can anticipate the applications of these properties by considering a 
/c-basis Bd C Id with distinct lead terms, as in the proof of Proposition 
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With respect to the lexicographic order, B,ir\k[xi, . . . , Xn] is then a A;-basis for 
Id n k[xi, . . . , Xn] for each i. With respect to the reverse lexicographic order, 
Bd n (xn) is then a /c-basis for Id n Thus, these orders enable us to find 
polynomials in an ideal which do not involve certain variables, or which are 
divisible by a certain variable. For a given degree d, one could construct such 
a basis Bd by applying Gaussian elimination to an arbitrary /c-basis for J^. 
However, this cannot be done for all d at once; such a computation would 
be infinite. We will finesse this difficulty by instead constructing a finite set 
of elements of I whose monomial multiples yield polynomials in I with every 
possible lead term. 

Such sets can be described as follows: 

Definition 2.4 A list F = [fi, . . . , fr] C / is a (minimal) Grobner basis for 
I i/in(/i), . . . ,in(/r) (minimally) generate in(/). 

in (7) is finitely generated because S is Noetherian, so Grobner bases exist 
for any ideal I. 

The order of the elements of F is immaterial to this definition, so F can 
be thought of as a set. We are using list notation for F because we are 
going to consider algorithms for which the order of the elements is signifi- 
cant. For convenience, we shall extend the notation of set intersections and 
containments to the fists F. 

A minimal set of generators for an ideal I need not form a Grobner 
basis for I. For example, if 5* = k[x,y] and / = + y'^,xy), then with 
respect to the lexicographic order, in(x^ + y"^) = x"^ and in{xy) = xy. Yet 
+ y'^) — x{xy) = y^ & I, so y^ G in(/). Thus, any Grobner basis for / 
must include y^; it can be shown that in(/) = (x^, xy, y^) and [x^ + y^, xy, y^] 
is a Grobner basis for I. 

On the other hand. 

Lemma 2.5 If F — [fi, . . . , f^] is a Grobner basis for I, then fi, . . . , fr 
generate I. 

Proof. For each degree d, we can construct a /c-basis Bd C Id with distinct 
lead terms, whose elements are monomial multiples of /i, . . . , Z^: For each 
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x"^ G in(/)rf, x"^ is a scalar multiple of x'^in(/j) for some x*-^ and some i; 
include fi in the set B^. Thus, the monomial multiples of fi, fr span 
/. I 



Proposition 2.6 (Spear, Trinks) Let R G S be the suhring R = 
k[xi, . . . ,Xn]- If F = [fi, . . . , fr] is a Grobner basis for the ideal I with 
respect to the lexicographic order, then FCiR is a Grobner basis for the ideal 
I n R. In particular, F (1 R generates I (1 R. 



Proof. (I PP77I , |Zac7§], Let / G / Hi?; in(/) is a multiple of in(/J 



for some i. Since in(/) G R, in(/j) G R, so ft G R. Thus, Ffli? is a Grobner 
basis for I n R. By Lemma F n R generates I H R. ■ 



Proposition |2.6| has the following geometric application: If / defines the 
subscheme X C P", then / fl k[xi, . . . , x„] defines the projection of X to 
P"-^ = Proj(A;[x„...,a:„]). 

Recall that the saturation P^^ of / is defined to be the largest ideal 
defining the same subscheme X C P" as I. P'"^^ can be obtained by taking 
an irredundant primary decomposition for /, and removing the primary ideal 
whose associated prime is the irrelevant ideal (xq, . . . / is saturated if 

J ^sat 

If the ideal I is saturated, and defines a finite set of points X C P", then 
/ n k[xn-i,Xn] is a principal ideal (/), where {/ = 0} is the image of the 
projection of X to P^ = Proj(/c[xn_i, Xn]). Given a linear factor of / of the 
form {bxn-i — aXn), we can make the substitution Xn-i = az, Xn = bz for 
a new variable z, to obtain from / an ideal J C k[xo, ■ ■ ■ , Xn-2, z\ defining a 
finite set of points in P"~^. For each point (cq, . . . , Cn-2, d) in the zero locus 
of J, (co, . . . , c„_2, ad, bd) is a point in the zero locus of /. 

If X C P"^^ is of dimension 1 or greater, then in general /n/i;[x„_i, x„] = 
(0), because a generic projection of X to P^ is surjective. In this case, an 
arbitrary substitution = bz can be made, and the process of 

projecting to P^ iterated. Thus, the lexicographic order can be used to find 
solutions to systems of polynomial equations. 

Recall that the ideal quotient (/ : /) is defined to be the ideal { G 
S\fg G /}. Since S is Noetherian, the ascending chain of ideals (/ : /) C 
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(/ : P) C (/ : /^) C ... is stationary; call this stationary limit (/ : f°°) = 
{g E S \ f^g G / for some m }. 

Proposition 2.7 // [x'^^ fi, . . . , x^^fr] is a Grobner basis for the ideal I with 
respect to the reverse lexicographic order, and if none 0//1, ■ ■ ■ , fr o,re divisible 
by Xn, then F = [fi, . . . , fr] is a Grobner basis for the ideal (/ : xjf ). In 
particular, /i, . . . , f,. generate (/ : x'^). 



Proof. ( ||Bay82|| , ||BS87a|| ) We have F C (/ : x^). Let / G (/ : x^); x'^f G / 



for some m, so in(x™/) is a multiple of in(x^'/j) for some i. Since ft is not 
divisible by Xn, in(/i) is not divisible by x„, so in(/) is a multiple of in(/j). 
Thus, F is a Grobner basis for (J : xjf ). By Lemma ^Tsj fi,...,fr generate 
(J : I 

If / = qo n qi n . . . n qt is a primary decomposition of /, then (/ : x^) = 
(riqi : xj^) = n(qj : xjf ). We have (q^ : xjf ) = (1) if the associated prime 
Pi of qi contains x„, and (q^ : xjf ) = q^ otherwise. Thus, if I defines the 
subscheme X C P", then (J : x^) defines the subscheme consisting of those 
primary components of X not supported on the hyperplane {x^ = 0}. 

(/ : xjf ) is saturated, because it cannot have (xq, . . . , x„) as an associated 
prime. If x„ belongs to none of the associated primes of / except (xq, . . . , x^,), 
or equivalently if {x„ = 0} is a generic hyperplane section of X C P", then 
(J : xjf ) = P^^. Thus, the reverse lexicographic order can be used to find the 
saturation of I. 

One of the most important uses of Grobner bases is that they lead to 
canonical representations of polynomials modulo an ideal /, i.e. a division 
algorithm in which every / G S* is written canonically as / = Oifi+h, where 
[/i, . . . , fr] is a Grobner basis for /, and h is the remainder after division. 

Recall the division algorithm for inhomogeneous, univariate polynomials 
/(x), g{x) G k[x]: Let in(/) denote the highest degree term of /. The 
remainder of g under division by / can be recursively defined by 

Rf{g) = Rf{g-cx'^f) 

if in(/) divides in(5f), where cx" = m{g)/m{f), and by 

Rfi9)=9 
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otherwise. 

Division can be generalized to homogeneous polynomials /i, ...,/,,, (7 G S, 
given a multiphcative order on 5" ( ||Hir64| | , | |Bri73| | , | |Gal74| | , | [Sch8U| |): The re- 
mainder Rpig) of g under division by the list of polynomials F = [fi,...,fr] 
can be recursively defined by 

Rp{g) = RF{g-c^^f,) 

for the least i so in{g) is a multiple cx"^ of in(/j), and by 

Rpig) = in{g) + Rpig - m{g)) 

if in(g) is not a multiple of any in(/j). Rpig) is an element of 5*. 

Thus, the fate of in{g) depends on whether or not in(5f) e 
(in(/i), . . . , in(/r)). Let / be the ideal generated by /i,...,/^. If F = 
[/i, . . . , fr] fails to be a Grobner basis for /, then the remainder is poorly 
behaved. For example, with respect to the lexicographic order on k[x,y], 

R[xy,x2+y2]{x'^y) = x^y - x{xy) = 0, 



but 



R[x^w,xy]i^'^y) = ^'^y - vi^'^ + y'^) 



so the remainder Rpig) is dependent on the order of the list F. Note that 
x'^y G (a;^ + y^, xy). 

If on the other hand, F is a Grobner basis for the ideal /, then Rpig) is a 
fc-linear combination of monomials not belonging to in(/). By Corollary p.3|, 
these monomials form a fc-basis for S/ 1, so each polynomial in S has a unique 
representation in terms of this fc-basis, modulo the ideal /. The remainder 
gives this unique representation, and is independent of the order of F (but 
dependent on the multiplicative order chosen for the monomials of S). In 
particular, Rpig) = if and only ii g E L 

An algorithm for computing a Grobner basis for / from a set of generators 
for / was first given by Buchberger ( [|Buc65|l , ||Buc76|| ). This algorithm was 



discovered independently by Spear ( [^pe77| 



, Bergman Per78|| , and 



Schreyer [|Sch80|| . It was termed the division algorithm by Schreyer, after the 
division theorem of Hironaka ( ||Hir64|| , ||Bri73|| , |pal74|| ) . 
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Define S{fi, fj) for i < j hj 

wliere x"^ = 6x^in(/j) = cx'-^in(/j) is tlie least common multiple of in(/j) 
and in(/j). 6x^/j and cx'"/j each have x^ as lead term, so x^ cancels out 
in S{fi, fj), and x^ > in(S'(/„ fj)). 

If F is a Grobner basis for the ideal J, then RF[S{fi, fj)) = for each 
i < j, since S{fi, fj) G /. Conversely, 



Proposition 2.8 (Buchberger) If Rp^S^fi, fj)) = for each i < j, then 
F = [fi, . . . , /r] is a Grobner basis for the ideal I = (/i, . . . , fr). 



See [p3uc65|] , ||Buc76|l . We postpone a proof until the theory has been 



extended to S-modules. This result can also be thought of as an explicit 
converse to the assertion that if F is a Grobner basis, then division is inde- 
pendent of the order of F: Whenever we have a choice in division between 
subtracting off a multiple of fi and a multiple of fj, the difference is a mul- 
tiple of S{fi,fj). If division is independent of the order of F, then these 



differences must have remainder zero, so by Proposition p.8| , F is a Grobner 
basis. 

As sketched in Section |l]. Proposition ^]8| can be used to compute a 
Grobner basis from a set of generators /i, . . . , for the ideal /: For each 
i < j so fr+i = RpiSifi, fj)) ^ 0, adjoin fr+i to the list F = [/i, . . . , /,.]. 
Note that fr+i G /. By iterating until no new polynomials are found, a 
Grobner basis F is obtained for /. This process terminates because S is 
Noetherian, and each new basis element corresponds to a monomial not in 
the ideal generated by the preceding lead terms. 

We now extend this theory to S'-modules. Let M be a graded, finitely 
generated S'-module, given by the exact sequence of graded S'-modules 



Ml Mo — > M — >0 



where Mq = Scqi © ... © Scog and Mi = Scu © ... © Scir are free S'-modules 
with deg(ejj) = dij for each i, j. We now think of F both as a list [/i, . . . , /r] 
of module elements, and as a map between free modules: Let fi = F{eii) ^ 
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for i = 1, . . . ,r, and let I C Mq be the homogeneous submodule generated 
by/i,...,/,. Thus, M = Mo//. 

A monomial of Mq is an element of the form x'^eoi; such an element has 
degree deg(x^) + doi. An order on the monomials of Mq is multiplicative if 
whenever x^^eoi > x^eoj, then x*-^x"^eoi > x'-^x^coj for all x'-^ G 5*. For some 
applications, such as developing a theory of Grobner bases over quotients of 
S, one wants this order to be compatible with an order on S: If x"^ > x-^, 
then one wants x'^eoi > x^eoi for i = 1, . . . ,r. The orders encountered in 
practice invariably satisfy this second condition, but it does not follow from 
the first, and we do not require it here. 

One way to extend a multiplicative order on S* to a compatible multi- 
plicative order on Mq is to declare x^eoi > x^eoj if i < j, or if z = j and 
-^A ^ -^B ^ Another way is to assign monomials x*^!, . . . , x'^'' in S to the basis 
elements Cqi, . . . , Cog of Mq, and to declare x^eoi > x^eoj if x^"*"*^' > x^+*-^^ , 
or if A + Ci = B + Cj and i < j. 

Fix a choice of a multiplicative order > on Mq. The constructions devel- 
oped for S carry over intact to Mq, with the same proofs ( ||Gal79|| , |[Sch80|] , 
[ Pay82|| ): Given an element / G Mq, define in(/) to be the lead term of 
/. Define in(/) to be the submodule generated by the lead terms of all el- 
ements of J C Mq] in(/) is a monomial submodule of Mq with the same 
Hilbert function as I. Define F = [fi, . . . , fr] C / to be a Grobner basis for 
I if in(/i), . . . ,in(/r) generate in(J); a set of generators for I need not be a 
Grobner basis for I, but a Grobner basis for / generates I. Given an element 
g G Mq, define Rpid) ^ Mq exactly as was done for the free module S. If F 
is a Grobner basis for /, then Rpig) = if and only if G /. 

The quotient of g under division by /i, . . . , can be recursively defined 

by 

Qrig) = cx'^eii + Qpig - cyi^fi) 
for the least i so '\n.[g) is a multiple cx"^ of in(/j), and by 

QF{g) = QF{g -Mg)) 

if in{g) is not a multiple of any in(/j). The quotient is an element of Mi. 

Following the recursive definitions of the remainder and quotient, it can 
be inductively verified that 

g = F{QF{g))+RF{g). 
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If F is a Grobner basis for /, and g & I, then Rpig) = 0, so the quotient 
hfts g to Ml. In this case, the quotient can be thought of as expressing g in 
terms of /i,...,/^. 

Define S{fi, fj) for i < j hy 

5(/„/,) = 6x^/,-cx^/„ 

if in(/j) and in(/j) have a least common multiple x^eok = &x^in(/j) = 
cx'^in(/j). Leave S{fi,fj) undefined if in(/j) and in(/j) lie in different sum- 
mands of Mq, and so don't have common multiples. 

Recall that the module of syzygies of /i, . . . , is defined to be the kernel 
of the map F, which is the submodule of Mi consisting of all h G Mi so 
F{h) = 0. Thus, ii h = hieii+. . .+hreir is a syzygy, then hifi+. . .+hrfr = 0. 
Let J C Ml denote the module of syzygies of fi,...,fr, and let K C Mi 
denote the module of syzygies of in(/i), . . . , in(/r). 

Define the map in(F) : Mi — ^ Mq by in(F)(eij) = in(/j); K is the kernel 
of in(F). For each i < j so S{fi, fj) is defined, define tij to be the element 

tij = fex^Cij — cx'^Cij G Ml, 

where x'^eofc = &x-^in(/j) = cx'-^in(/j) is the least common multiple of in(/j) 
and in(/j), as before. m{F){tij) = 0, so each t^j belongs to the syzygy module 
K. Observe that F{tij) = S{fi, fj). 

Assign the following multiplicative order on Mi, starting from the order 
on Mo ([gcESg; see also ||MM86|| ): Let x^di > x^e^ if x^in(/i) > x^in(/j). 



or if these terms are /c- multiples of each other and i < j. If the order on Mq 
is compatible with an order on S, then this order on Mi is compatible with 
the same order on S. 

With respect to this order on Mi, we have 

Lemma 2.9 The list [tij] is a Grobner basis for the module K of syzygies 
o/in(/i),...,in(/r). 

Proof. Let h e Mi, so 'm[F)(h) = 0. Then in(F)(in(/t)) is canceled by 
in(F)(/i — in(/i)) in Mq. Therefore, if m{h) = x'^eij, then h has another term 
x^Cij so x^in(/j) and x^in(/j) are /c-multiples of each other and i < j. 
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Thus, tij is defined and in(tjj) divides in(/i), so [tij] is a Grobner basis for 
K. I 

Tlius, tlie set {tij} generates K. In general, tlie [Uj] are far from being 
a minimal Grobner basis for K; we consider the effects of trimming this list 
in Proposition |2.1U| below. 



Define 

Sij = tij — QpiS^fi, fj)) 

whenever Rp^S^fi, fj)) = 0. Note that in(sjj) = in(tij). Each Sij is the 
difference of two distinct elements of Mi, each of which is mapped by F to 
S{fi, fj), so F{sij) = 0. In other words, Sij belongs to the syzygy module J. 
Conversely, 



Proposition 2.10 (Richman, Spear, Schreyer) Choose a set of pairs 
T = {{i,j)} such that the set {tij}{i,j)eT generates the module K of syzy- 
gies o/in(/i), . . . ,in(/^). If RpiSifi, fj)) = for each {i,j) G T, then 

(a) F = [fi, . . . , fr] is a Grobner basis for I ; 

(b) the set {sij}(ij)gr generates the module J of syzygies of fi, . . . ,fr- 
Moreover, 

(c) ^/ [^ij ](i,i)er is a Grobner basis for K , then [sjj](jj)gT is a Grobner 
basis for J. 



Proof. ( lEicTl , |Spe77|] , ||Zac78|] , [ gch80| ) First, suppose that [tij](i,i)eT is 
a Grobner basis for K. Let h & J, so F{h) = 0. By the same reasoning as in 
the proof of Lemma |2.9| , we can find {i,j) G T so in(tjj) divides m{h). Since 
in(sjj) = in(tjj), in(sjj) also divides in(/i), so [sij](ij)^T is a Grobner basis 
for J, proving (c). 

Now, suppose that {tij} {i,j)<^T merely generates K. Let T' be a set of 
pairs so [t^m ](£,m)GT' is a Grobner basis for K. It is enough to construct a list 
[uim](£,m)&T' of elements of J, generated by {sij}(ij)(zT, so m{uim) = m{tim) 
for all (£, m) G T'. Then by the preceding argument, [uim]{e,m)<=T' is a 
Grobner basis for J, so {sjj}(jj)gT generates J. 

Write each tim = J^gemijtij, for (£, m) G T' and (i, j) G T, in such a way 



that the terms of tem and each term of each product gimijiij map via in(F) 
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to multiples of the same monomial in Mq. In other words, find a minimal 
expression for each tem, which avoids unnecessary cancellation. Then define 



We have in(-u^m) = in(t^m), proving (b). 

Let f & I, and choose g G Mi so / = F{g). Let h G Mi be the remainder 
of g under division by [uim](£,m)eT'] f = F{h). Since m{h) is not a multiple 
of any m{uem) = the lead term of F{m{h)) is not canceled by any 

term of F{h — in{h)). Therefore, if m{h) = ax^eij, then in(/j) divides in(F). 
Thus, F = [fi, . . . , fr] is a Grobner basis for /, proving (a). ■ 



Proposition p.8| follows as a special case of this result. 

The above proof can be understood in terms of an intermediate initial 
form ino(/i) for h G Mi: Apply the map in(F) separately to each term of h, 
and let x"^ G Mq be the greatest monomial that occurs in the set of image 
terms. Define ino(/?.) to be the sum of all terms of h which map via in(F) 
to multiples of x"^. Then in refines Iuq, for according to the order we have 
defined on Mi, m{h) is the term of ino(/?.) lying in the summand of Mi whose 
basis element Cj has the smallest index i. 

In this language, tij = ino(tij) = ino(sjj). Our expressions for the tirn 
have the property that each ggmijtij = i^o{gemijtij), with each term of each 
product for a given tgrn mapping via in(F) to multiples of the same monomial 
x^. Thus, each ino{gemijSij) = gemijUj] the tails gemijisij - ino(sij)) stay out 
of our way, mapping termwise via in(F) to monomials which are less than 
x"^ with respect to the order on Mq. 

Observe that Qrig) is a linear combination of monomials not belonging 
to in(J), for any g G Mq. 

In ||Buc79|| , Buchberger gives a criterion for selecting a set T of pairs 



in the case where I is an ideal: If (io, ^i), (h, ^2), ■ ■ ■ , (^s-i, is) £ T, and 
the least common multiple of in(/j(,), in(/jj, . . . , in(/jj is equal to the least 
common multiple of iia{fi^) and in(/jj, then {iq, is) need not belong to T. In 
other words, if ti^i^ G (tjoii, . • . then the pair {io,is) is unnecessary; 

this condition is equivalent to the condition of Proposition p.lO| , for the case 
of an ideal. 

Suppose that we wish to compute the syzygies of a given set of elements 
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gi, . . . ,gs of Mq. To do this, compute a Grobner basis /i, . . . , for the 
submodule I C Mq generated hj gi, . . . , gs- Keep track of how to write each 
fi in terms of gi, . . . , gs- Using these expressions, each syzygy of /i, . . . , 
can be mapped to a syzygy of gi, . . . , gg- These images generate the module 
of syzygies of gi, ■ ■ ■ ,gs', the set of syzygies obtained in this way is not in 
general minimal. 

Syzygies can be used to find a minimal set of generators for a submodule 
/ C Mq from a given set of generators gi, . . . , gs'. If higi + . . . + hrgr = is 
a syzygy of ^fi, . . . , gs with hi e k, then gi = (/i2fi'2 + • • • + Kg,) /hi, so gi is 
not needed to generate /. All unnecessary generators can be removed in this 
way. 

Alternatively, a careful implementation of Grobner bases can directly find 
minimal sets of generators for submodules: Starting from an arbitrary set 
of generators, we can eliminate unnecessary generators degree by degree, by 
removing those which reduce to zero under division by a Grobner basis for 
the ideal generated by the preceding generators. 

Either way, we can trim the set of syzygies computed via Grobner bases for 
a given set of generators (71, . . . , of J, to obtain a minimal set of generators 
for the syzygy module J. By starting with a minimal generating set for /, 
and iterating this method, a minimal free resolution can be found for I. 

A beautiful application of these ideas yields a proof of the Hilbert syzygy 
theorem, that minimal free resolutions terminate (Schreyer ||Sch80|| , [^chQl 



for an exposition see also Eisenbud ||Eis92|| ). At each stage of a resolution. 



order the Grobner basis F for / in such a way that for each i < j, letting 
in(/j) = ax'^eofc and in(/j) = 6x^eo^, we have x"^ > x^ in the lexicographic 
order. If the variables xi, . . . , Xm are missing from the initial terms of the fi, 
then the variables xi, . . . , Xm+i will be missing from the initial terms of the 
syzygies Sij. Iterating, we run out of variables, so the resolution terminates. 



3 Bounds 

How hard are the algorithms in algebraic geometry? We describe some key 
bounds. The best known example is the bound established by G. Hermann 



Her26| for ideal membership: 
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Theorem 3.1 (G. Hermann) Let k be any field, let {fi,...,fk) C 
k[xi, ...,Xn] and let d = max(deg(/i)). If g & (/i,...,/^), then there is an 
expression 

k 

9 = Y1 ^ifi 



where deg(aj) < deg{g) + 2{kd) 



i=l 



This type of bound is caUed "doubly exponential" . However, with the ad- 
vent of the concept of coherent sheaf cohomology |[Ser55|| and the systematic 



study of vanishing theorems, it has become apparent that the vanishing of 
these groups in high degrees is almost always the most fundamental bound. 
The concept of an ideal being "m-regular" or "regular in degrees > m" was 
introduced by one of us ||Mum66|| by generalizing ideas of Castelnuovo: 



Definition 3.2 Q Let k be any field, let I C k[xo, . . . be an ideal gener- 
ated by homogeneous polyomials, let Id be the homogeneous elements in I of 
degree d, let I be the corresponding sheaf of ideals in Opn, and let I{d) be 
the d*'^ twist of I. Then the following properties are equivalent and define the 
term "m-regular" : 

(a) the natural map Im ~^ II^{I{m)) is an isomorphism and W{I{m — i)) 
= (0), l<i<n 

(b) the natural maps Id II^{X{d)) are isomorphisms for all d> m and 
H'{I{d)) ={0)ifd + i>m,i>l. 

(c) Take a minimal resolution of I by free graded k[X]-modules: 
0^ © kM-Can ^ ...^ © A;[x] -e^o ^ A;[x] ^ A;[x]//^0. 

Then deg(ea,j) < m -\- i for all a, i. (In particular, if fa = (t>o{ea,o), then 
/i, . . . , frg are minimal generators of I, and deg(eQ,_o) = deg(/Q) < m.) 

The intuitive idea is that past degree m, nothing tricky happens in the 
ideal /. Unfortunately, neither (a), (b) nor (c) can be verified by any obvious 
finite algorithm. This lack of a finitely verifiable criterion for m-regularity has 
been remedied by a joint result of the first author and M. Stillman ||BS87a| ]: 

^ The definition has been slightly modified so as to apply to ideals / instead of the 
corresponding sheaf of ideals X. 
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Theorem 3.3 (Bayer-Stillman) / is m-regular if and only if the degrees 
of the minimal set of generators of I are at most m, and there exists a set 
Ho, . . . ,y£ of linear combinations of Xq, ...,Xn such that for all homogeneous f 
of degree m, 

yofel =^ fel 

yif el ^ f el + /c[x] • yo 

e-i 

ytfel f e I + J2k[x]-yi 

i=0 

and ^ 

feI + Y.k[^]-yi. 

i=0 

Moreover, if this holds at all, it holds for yo, . . . ,yi taken arbitrarily from a 
Zariski-open set in the space of £ + 1 linear forms. 

To see why m-regularity is a key bound, we want to show that it controls 
some of the geometric features of the ideal /. Let's introduce several refined 
notions of the "degree" of /: 

Definition 3.4 //" 7 = qo n qi fl . . . fl Z5 a primary decomposition of 
I J V^Qi — Pi prime and V"(p,;) is the suhvariety Zi ofP'^ for i > 1, while 
Po = (^^O) • • • I ^n) (so that V{po) = then first let qi, qs be the isolated 
components, (i.e., Zi Zj if 1 < i < s, 1 < j < t, i ^ j, or equivalently, 
V{I) = Zi U . . . U Zs is set-theoretically the minimal decomposition of V{I) 
into varieties). Then let 

mult(qj) = length £ of a maximal chain of Pi-primary ideals: 
qi = Je'^ Ji-i ^ ■ ■ ■ ^ Ji — Pi 

(Equivalently, this is the length of the local ring A;[x]p./7A;[x]p., or, in the 
language of schemes, if r] is the generic point of Zi, then this is the length of 

O pnJ 

deg(Zj) = usual geometric degree of Zi'. 
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the cardinality of Zif] L for almost all 
linear spaces L of complementary dimension. 



geom-deg^(/) 



mult(qi)deg(Zj) 



i such that dim Zi = r 



l<i<s 



If Qj is one of the non-isolated, or embedded components, then we extend 
the concept of multiphcity more carefully: Let 



li = |nqj \j such that Pj ^ or equivalently Zj ^ Z^j n 
and 

mult7(qj) = length £ of a maximal chain of ideals: 
n /j = ^ Je-i ^ . . . ^ Jo = /j 
where each satisfies: ab E J^, a ^ p, =^ 6 e J^. 

(Equivalently, Jk equals qfcn/j for some pj-primary ideal q^.) In particular: 



For s-|-l<i<t, an equivalent way to define mult7(qj) is as the length of 
the module 




mult7(qo) 








or, in the language of schemes, the length of 



Pi 



where rj is the generic point of Z^. 
Then write 



arith-deg^(/) 



E 



mult/(qi) deg(Zj) 



i such that dim Zi = r 
1 < i < s 
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and 

arith-deg_]^(/) = mult/(qo). 
The idea here is best illustrated by an example: let 

I = {x\,XiX2) C k[xQ,Xi,X2]. 



Then 



Then 



so 



I = qi n q2 

qi = [x]), pi = (xi), Zi = {line xi = 0} 

?2 = (xi,xiX2,X2), P2 = (a^^i, 2^2), -^2 = {point (1,0,0)}. 

deg(Zi) = 1, mult(gi) = 1 



geom-deg]^(/) = arith-deg]^(/) = 1. 
One might be tempted to simply define 
mult/(q2) = length of chain of p2-primary ideals between q2, P2 
and since 

A;[x]q2/q2A;[x]p2 ^ K ■ I + K ■ xi + K ■ X2,K ^ k{xo) 
this is 3. But embedded components are not unique! In fact, 

I = qinq'2 

Q2 — (^^1^2) also, 

which leads to 

k[^]pJq'^k[^]p,^K-l + K-X2 

which has length 2. The canonical object is not the local ring 
/c[x]p2/q2/i:[x]p2 but the ideal 

Ker {k[x.]pjlk[x.]p^ A;[x]p2/p2/c[x]pJ ^ k ■ xi 
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which has length 1. Thus, the correct numbers are 



muh/(q2) = 1 

and 

geom-degQ(/) = 
arith-degQ(/) = 1. 

Now the question arises: find bounds on these degrees in terms of gen- 
erators of I. For geometric degrees, a straightforward extension of Bezout's 
theorem gives: 

Proposition 3.5 Let d{I) be the maximum of the degrees of a minimal set 
of generators of I . Then 

geom-deg,(/) < rf(J)"-^ 



A proof can be found in ||MW83|| . The idea is clear from a simple case: 



Suppose f,g,h G K[x, y, z] and f = g = h = consists of a curve C and £ 
points Pi off C. We can bound i like this: Choose 2 generic combinations /', 
g' of /, g, h so that f' = g' = does not contain a surface. It must be of the 
form C U C, C one-dimensional, containing all the Pi but not the generic 
point of C . Then by the usual Bezout theorem 

degC'<deg/' degg' = d{lf. 

Let h' be a S^d generic combination of /, g, h. Then C fl {h' = 0} consists of 
a finite set of points including the Pj's. Thus 

i = m<i^iC'n{h' = o}) 

< degC' ■ d{I) by Bezout's theorem 

< d{I)\ 

Can arith-deg(/) be bounded in the same way? In fact, it cannot, as we will 
show below. Instead, we have 
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Proposition 3.6 If m{I) is the regularity of I , then for — 1 < r < n, 

arith-deg,(/) < + n - r - l\ 

\ n — r J 

which replaces d{I) by the regularity of /. A proof is given in the technical 
appendix. 

We have introduced two measures of the complexity of a homogeneous 
ideal I. The first is d{I), the maximum degree of a polynomial in a minimum 
set of generators of /. The second is m{I), which bounds the degrees of gen- 
erators and of all higher order syzygies in the resolution of / (Definition 
(c)). Obviously, 

d{I) < m{I). 

A very important question is how much bigger can m(/) be than d{I)l 
The nature of the answer was conjectured by one of us in his thesis Pay82|| 
and this conjecture is being borne out by subsequent investigations. This 
conjecture is that in the worst case m(/) is roughly the (2")*^ power of d{I) 
- a bound like G. Hermann's. But that if / = I{Z) where Z is geometrically 
nice, e.g. is a smooth irreducible variety, then m(/) is much smaller, like the 
n^J^ power of d{I) or better. This conjecture then has three aspects: 

(1) a doubly exponential bound for m(J) in terms of d{I), 
which is always valid, 

(2) examples of / where the bound in (1) is best possible, or nearly so, 

(3) much better bounds for m(/) 

valid if V{I) satisfies various conditions. 

All three aspects are partially proven, but none are completely clarified 
yet. We will take them up one at a time. 

A doubly exponential bound for m(J) in terms of d{I) may be deduced 
easily in characteristic zero from the work of M. Giusti ||Giu84|| and A. Galligo 

Theorem 3.7 //char(A;) = and I C fc[xo, . . . is any homogeneous ideal, 
then 

m{I) < {2d{I)f"'\ 
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It seems likely that Theorem holds in characteristic p, too. A weaker 
result can be derived quickly in any characteristic by straightforward coho- 
mological methods: 

Proposition 3.8 If I C k[x(), . . . x„] is any homogeneous ideal, then 

m{I) < {2d{I) f-. 

The proof is given in the technical appendix. 

Next, we ask whether Theorem |3.7| is the best possible, or nearly so. The 
answer is yes, because of a very remarkable example due to E. Mayr and A. 
Meyer ||MM82|| . 



Example 3.9 Let be the ideal in lOn variables S^"'\ F^™), cf™\ 1 < 
i < 4,1 < m < n defined by the lOn — 6 generators 

2<m<n i - S^^-i^Cj™"^) 



p{m)(j{-nT-) _ p{m)(j{m) 
1 < m < n — I < g(m) (jijri) j^{m.) _ g{m)(j(m) j^M 

_ pi^ri) j^{rn) 

Let be the ideal gotten from by homogenizing with an extra variable 
u. Then Mayr and Meyer ||MM82|, lemma 8, p. 318] prove: 



Lemma 3.10 Let e„ = 2^". // M is any monomial in these variables, 

and ^ - 5(")M G if and only if 

M = \ 
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Now note that the generators of and are all of the very simple type 
given by a difference of two monomials. Quite generally, if 

J C • • • , Xn\ 



J = (...,x"* -x^%...) 



l<i<fe 



then the quotient ring A;[x]/J has a very simple form. In fact, we get an 
equivalence relation between monomials generated by 

x'^i+T x^*+^, any i,-f 

and 

A;[x]/J = ©5 A; ■ x^ 
where 6 runs over a set of representatives of each equivalence class. 

Bearing this in mind, let's look at the 1^^ order syzygies for the homoge- 
neous ideal: 

5"^") and F^"^ are part of a minimal set of generators, and let G 
complete them. Then syzygies are equations: 

One such is given by: 



for some and some e > (the extra power is necessary because some 
terms Rafa have degree greater than e„ + 2) whose degree is 2 + e„ + e. Now 
express this syzygy as a combination of a minimal set of syzygies. This gives 
us in particular: 

^e„+e^(n) ^ ^^^^^ 

Then for some A, px must have a term of the form or C'i"'\ hence 
the monomial UiS^"'^ or Ui 

^(n) g(n) 

occurs in pxS^"'\ But by the general 
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remark on quotient rings by such simple ideals, this means that this term 
must equal some second term M S'*-"-' (M a monomial in p\) or M F^'^^ (M 
a monomial in q\) mod . By the lemma, the first doesn't happen and the 
second only happens if the term C\^^ {bI^'^Y" occurs in qx, in which case 
e„ + 1 < deg qx = deg(syzygy(pA, q\, Rax)) - 1- This proves: 

Proposition 3.11 has for its bounds: 

d{J) = 4 
m(J) > 2'" + l. 

Going on to the 3^*^ aspect of the conjecture, consider results giving better 
bounds for m{I) under restrictive hypotheses on V{I). 

Theorem 3.12 If Z G P" is a reduced subscheme purely of dimension r, 
and I = I{Z) is the full ideal of functions vanishing on Z , then 

(a) if r < 1, or Z is smooth, char{k) = and r < 3, then: 

Tn{I) < deg Z — n + r + 1 

(b) if char{k) = and Z is smooth, 

m{I) < (r + l)(deg(Z) - 2) + 2. 



Since deg(Z) < (Proposition |3.5|), these bound m(/) in terms of 



d{I). 



1, and to Pinkham 

3. It is conjectured by Eisenbud and Goto |[EG84 



Part (a) of this are due to Gruson-Lazarsfeld-Peskine ||GLP83|| for r < 
Pin86 |, Lazarsfeld [ Laz87 |, and Ran [ Ran9C ] for r < 

and others, that the 



bound in (a) holds for all reduced irreducible Z, and it might well hold even 
for reduced equidimensional Z which are connected in codimension 1. As 
this problem is now understood, the needed cohomological arguments follow 
formally, once one can control the singularities of a projection of the variety. 
These singularities become progressively harder to subdue as the dimension 
of the variety increases, and are what impedes definitive progress beyond 
dimension 3. 
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Part (b) is due to the second author and is proven in the technical ap- 



pendix. It has been generahzed by Bertram, Ein, and Lazarsfeld ||BEL91 
to show that any smooth characteristic variety of codimension e de- 
fined as a subscheme of P" by hypersurfaces of degrees di > . . . > dm is 
{di + . . . de — e + l)-regular. Since we cannot decide the previous conjecture, 
this is a resuh of considerable practical importance, for it strongly bounds the 
complexity of computing Grobner bases of smooth characterisitic varieties 
in terms of the degrees of the input equations. 

The biggest missing link in this story is a decent bound on m(J) for any 
reduced equidimensional ideal /. We would conjecture that if a linear bound 
as in part (a) doesn't hold, at the least a so-called "single exponential" bound. 
I.e. m (/) < ci°(") ought to hold. This is an essential ingredient in analyzing 
the worst-case behavior of all algorithms based on Grobner bases, and would 
complete the story about what causes the bad examples discussed above. At 



least in some cases Ravi | Rav90 ] has proven that the regularity of the radical 



of a scheme is no greater than the regularity of the scheme itself. 

There is a direct link between the bounds that we have given so far and 
the G. Hermann bound with which we started the section. This results from 
the following: 

Proposition 3.13 Let C k[xi, . . . , x„] have generators fi,---,fk and let 

C k[x(), Xi, . . . , Xn] be the ideal generated by homogenizations . . . , 
of the fi. Let = qo H . . . H qt be the primary decomposition of , let 
Zi = V{qi) and let 

multoo(/^) = max [mult/(qjj + . . . + mult/(qjj + mult/(qo)] 

where the max is taken over chains V{{xq)) 'D Zi^^ . . . ^ Zi^. If g & I^, then 
we can write: 

k 

9 = ^(^ifi 

4 = 1 

where 

degOi < degfif + multoo(/^). 

The proof goes like this: Let g^ be the homogenization of g. Consider the 
least integer m such that x'^g^ G . Since g & I, this m is finite. Moreover, 
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if 

then 
and 

degOi = deg(a^) < deg{x'^'g^) - deg fj < m + deg{g). 
Now in the primary decomposition of , suppose that for some k, 

x^g^ e n qi, and x^^^g^ ^ if j ^ S. 

Choose i ^ S such that V{qe) is maximal. Since g e I^, we know V{qe) C 
V{{xo)), hence xq G p^. Let 

Then mult/(q^) is easily seen to be the length of a maximal chain of ideals 
between: 

/ • k[x]p^ and Is ■ k[x]p^. 
But look at the ideals Jp, for p > 0, defined by 



If Jp — Jp+i, then 



Jp 



i.e., 4+V = a4+^+y + 6, 6e7. 
But 1 — axo is a unit in A;[x]p^, so 

Jp = = (1 - axo)-^b e /A;[x]p,. 

This means that in any case 
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hence, because qi is prpriniary: 



Induction now shows that 



fc+multj(q^) 

^0 9 ^ 



mnltooil") h ^ jH 



Corollary 3.14 Let I^, he as above. If g E I^, then 
where deg(ai) < deg{g) + 



m(I)+n+l 
n+1 



Proof. Combine Propositions p.6| and |3.13 . 



If we further estimate m{I) by Theorem p.7| in characteristic or by 
Proposition |3]^, we get somewhat weaker versions of Hermann's Theorem 
But if I = V (Z) , Z a good variety, we may expect the Corollary to give much 
better bounds than Theorem 

Corollary |3.14| shows that any example which demonstrates the necessity 
of double exponential growth in Hermann's ideal membership bound (The- 
orem |3.1| ) also demonstrates the necessity of double exponential growth in 
the bounds on m{I) given in Theorem |3.7| and Proposition p.8| . Thus we can 
make use of the general arguments for the existence of such examples given in 
[[MM82|| , rather than depending on the single example of Proposition |3.11| , to 



show that the bounds on m{I) inevitably grow double exponentially: Since 
in Corollary |3.14| , the degrees of the Oj are bounded by a single exponential 
function of m{I), in all examples where the degrees of the Oj grow double 
exponentially, m{I) also grows double exponentially. 

This line of argument gives a geometric link between the ideal membership 
problem and m(J): In Corollary p.l4| , if exhibits Oj of high degree, then 
has primary components of high multiplicity. These components force m(J) 
to be large, and distinguish from good ideals considered in Theorem |3.12| 
and related conjectures. 

A major step in understanding the gap between the double exponential 
examples and the strong linear bounds on the regularity of many smooth 
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varieties was taken by Brownawell [[Bro87|| and Kollar ||Kol88|| . They discov- 
ered the beautiful and satisfying fact that if we replace membership in / by 
membership in \fl^ then there are single exponential bounds on the degrees 
of a~: 



Theorem 3.15 (Brownawell, Kollar) Let k he any field, let I = 
(/i, •••,/fe) C k[xi, ...,Xn] and let d = max(deg(/i), i = 1, ■ ■ ■ , A;; 3). Ifn = l, 
replace d by 2d — 1. If g & Vl, then there is an expression 

k 
i=l 

where s < d^ and deg(aj) < (1 + dQg{g))d'^ . In particular: 



What this shows is that although the bad examples have to have primary 
components at infinity of high degree, nonetheless these primary ideals con- 
tain relatively small powers of vT^. The picture you should have is that 
these embedded components at infinity are like strands of ivy that creep a 
long way out from the hyperplane at infinity, but only by clinging rather 
closely to the affine components. 



Technical Appendix to Section |^ 



1. Proof of the equivalence of the conditions in Definition p.2| : 

In [|Mum66| , pp. 99-101], it is proven that for any coherent sheaf on 
P", H\T{-i)) = (0), z > 1 implies that the same holds for T{d), all ci > 0, 
and that H'^{T{d)) is generated by H'^{T) ® H^{0{d)). In particular, if you 
apply this to T = I{m), the equivalence of (a) and (b) follows. (Note the 
diagram: 

n n 

k[^]d — ^ H%Opn{d)) 

which shows that H^{I{m)) is injective for every d). To show that (b) 



(c), first note that we may rephrase the reults in ||Mum66|| to say that if 
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H^{J^{—i)) = (0), i > 1, then the degrees of the minimal generators of the 
/c[x]-module 

are all zero or less. So we may construct the resolution in (c) inductively: at 
the fcth stage, say 



a=l 



^[x] ■ e„,fc„i — > ■ ■ ■ — > fc[x] — > fc[x]/J — > 



has been constructed, let Mk = kei^cpk-e) and let J-'k be the corresponding 
sheaf of ideals. The induction hypothesis will say that H^{J-'k{m + — 1)) = 
(0), i > 1. Therefore Mk is generated by elements of degree < m + k, i.e., 
da = deg Cq fe < m + /c, all a. We get an exact sequence 



hence 



Mk+i © (A;[x] ■ e«,fc) — . il/4 ^ 



^ ^ © Opn(-d„) ^ Tk ^ Q (1) 



Therefore 



© HHOjinim+k-i-da)) — ^ HHj^kim + k-i)) — y (2) 

q:=1 

H'+\j^k+i{jn + {k + l)-{i + l))) — > © /7*+^(Opn(m+A;-i-d„)) 

a=l 

is exact. But m + A; — z — do, > — i so iJ*+-'^((9pn(m + A; — z — do,)) = (0). This 
shows that Tk+i satisfies the induction hypothesis and we can continue. Thus 
(c) holds. To see that (c) ^ (a), we just use the same exact sequences (1) 
and prove now by descending induction on k that Wi^J^kijn + k — i)) = (0), 
i > 1. Since / = jFg, this does it. The inductive step again uses (2), since 
H'(Opn{m + k-i- da)) = (0) too. 



2. Proof of Proposition |3.6| : 

Look first at the case r = 0. Let X be the sheaf of ideals defined by / 
and let I* D X be the sheaf defined by omitting all 0-dimensional primary 
components of /. Consider the exact sequence: 

— > X(m-l) — ^X*(m-1) — > (X7X)(m-l) — > 
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This gives us: 

/7°(X*(m-l)) — . H%{T*/I){m-l)) — > H\T{m - I)) 

Now H^{I{m — 1)) = (0) by m-reg ularity, and h°{{I*/I){m - 1)) = 
/i°(X*/X) = length(X*/X) = arith-dego(/) since X*/X has 0- dimensional sup- 
port. But H^il*{m - 1)) C //°(Opn(m - 1)), so 



If r > 0, we can prove the Proposition by induction on r. Let H be 
a generic hyperplance in P", given by /i = 0. Let Ih = {I,h)/{h) C 
k[xo, . . . , Xn]/{h) = k[x'Q, . . . , x'j^_i] for suitable linear combinations of Xi. 
Then it is easy to check that: 



arith-dego(/) < h° 



(X*(m-1)) 
(Cpn(m- 1)) 




arith-degr(/) = arith-degr„i(/H) 



and that Ih is also m-regular, so by induction 





m + n — r 



n — r 



1 



If r = — 1, we use the fact that 








if (i > m, hence 





3. Proof of Proposition |3.8|: 



39 



Let I G k[xQ, . . . Xn] and assume, after a linear change of coordinates, that 
Xn is not contained in any associated prime ideals of /. Let / C k[xQ, . . . Xn-i] 
be the image of /. Then d{I) — d{I) and by induction we may assume 

m(7) < (2(i(/))("-^)'. 

We will prove, in fact, that 

/ TV /-FX / m(I) — 1 + n\ 
m{I) < m{I) + ( ^ 1 (3) 

and then we will be done by virtue of the elementary estimate: 

if m* = (2ci(/))("-^)', and d > 2, then m* + ~J ^ j < (2d(/))"' 

To prove (3), we use the long exact sequence 
{I:(xo))k-i ^ h 7k 

H\X{k-l)) H'^mk-l)) H^'ilik-l)) ^ 

H\l{k-1)) H\X{k)) H\T{k)) 

where (7 : (xq)) ^ {f\xof e I}. Let m = m(7). Note that W{T{k - 1)) = 
(0), i > 1, A; > m, hence 

H'{I{k - 1)) ^ H\I{k)) 

is an isomorphism if/c>7n — 1 + 1 and i >2. Since H^{I{k)) = (0), A; 3> 0, 
this shows that H\l{k)) = (0), i>2,k>m-i. Moreover 7^ H^{X{k)) 
is an isomorphism if /c > m, hence 5 = ii k > Wi, hence H^{X{k)) = (0), 
k > fn — 1. But now look at the surjectivity of Ik — > H^{X{k)). For all 
k, let Mfe be the cokernel. Then ©M^ is a A;[x]-module of finite dimension. 
Multiplication by xq induces a sequence: 

^ _^ {I : {xo))k-i _^ ^ Mk ^ 
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which is exact if A; > m. But if, for one value of A; > m, 

(/ : {xo))k = h (4) 

then by Theorem |3.3| , / is fc-regular and (4) continues to hold for larger k, 
and Mk must be (0). In other words, 

dim Mk , k > m — 1 

is non increasing and monotone decreasing to zero when k>m. Therefore 

m{I) < 
< 

< 



m + dim M-, 



m—1 

TFT + dim A;[x]m_i 
— 1 + 
K n I 



m + 



which proves (3). 



4. Proof of Theorem 3.12| (b): 



Let Z he a smooth r-dimensional subvariety of P" and d = degree of Z. 
We first consider linear projections of Z to P** and to P''"^. To get there, 
let Li C P" be a linear subspace of dimension n — r — 1 disjoint from Z and 
L2 C Li a linear subspace of dimension n — r — 2. Take these as centers of 
projection: 

P"-Li D Z ^ 

pr+l _ |p| 3 2^ 
i 

P2 ^ pr 

Let XQ,...Xr+i be coordinates on P^+i so that p = (0, ...,0, 1), hence 
xq, • • • ^r t^re coordinates on P^. Let /(xq, . . . x^+i) = be the equation 
of the hypersurface Zi. 

Now there are two ways of getting r-forms on Z: by pullback of r-forms 
on P'' and by residues of (r + l)-forms on P**"^ with simple poles along Zi. 
The first gives us a sheaf map 
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whose image is -Bi), Bi the branch locus of p2. Corresponding to this 
on divisor classes: 

Kz = P*2{Kpr)+B^ (5) 

= -{r + l)H + Bi, 
where H = hyperplane divisor class on Z. The second is defined by 
, , dxi A ... A dxr+i * f / \ dxi A ... A dxr \ 

-f " i"''') ■ a//ax,« j 

and it gives us an isomorphism 

B2 is a divisor which can be interpreted as the conductor of the affine rings 
of Z over those of Zi: i.e., 

feOzi-B2) ^ f ■ iPi,.Oz) c Oz,. 

In particular, 

Pi,,(Oz(-52)) ^ sheaf of Oz, - ideals C in C^,. (7) 

A classical reference for these basic facts is Zariski ||Zar69|| , Prop. 12.13 
and Theorem 15.3. A modern reference is Lipman [[Lip84|| (apply Def. (2.1)b 
to pi and apply Cor. (13.6) to Zi C P''"'"^). (4) gives us the divisor class 
identity: 

KZ + B2 = pl{Kpr+l+Z,) (8) 

= {d-r-2)H. 

(5) and (8) together tell us that 

Bi + B2 = {d- l)H. 

In fact, the explicit description (6) of the residue tells us more: namely that 
if 2/1, . . . , l/r are local coordinates on Z, then 

A? ^ ■ aTTa ayi A . . . A dyr 
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generates ^^^(fi2) locally. But = is a local equation for Bi, so this 

means that df /dxr+i = is a local equation for Bi + B2. But df /dxr+i = 
is a global hypersurface of degree d — 1 in P^+^, hence globally: 

b, + b,^pUv(^)) 

(equality of divisors, not merely divisor classes). All this is standard classical 
material. 

(7) has an important cohomological consequence: let C* C Opr+i be the 
sheaf of ideals consisting of functions whose restriction to Zi lies in C. Then 
we get an exact sequence: 







pr+l 



-Zi] 



c* o 



pr+l 



CO 







hence an exact sequence 

^ Opr+i(e-d) ^ c*Oj,r+i{£) ^ pi,,{Oz(m - B2)) 

for all integers £. But H\Opr+i{i - d)) = (0), hence 

H%C*Opr+i{i)) ^ H\Oz{£H-B2)) 
is surjective, hence 

H\Oz{m-B2)) C Im[i7°(Op.+i(£)) ^ H\Oz{iH)) 



(9) 



Now let us vary the projections pi and p2. For each choice of Li, we get 
a different Bi. call it Bi{Li), and for each choice of L2, as different B2. call 
it B2{L2). By (5) and (8), all divisors Bi{Li) are linearly equivalent as are 
all divisors ^2(1/2). Moreover: 



n 5i(Li 







n B2{L2) = 

This is because, ii x & then there is a choice of Li such that pi : Z — > P'' 
is unramified at y] and a choice of L2 such that p2{x) £ Zi is smooth, hence 
P2 is an isomorphism near x. Thus 



Bi(Li)| = \Kz + {r + l)H\ 
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and 

\B2{L2)\ = \ Kz + id-r-2)H\ 
are base point free linear systems. 

Next choose (r + 1) L2's, called L2,l<a<r + 1, so that if i?2"' = 
B2(L'f^), then f] i?2"^ = 0- Look at the Koszul complex: 

a 

This is exact and diagram chasing gives the conclusion: 
W{Oz{iH + 1)B2)) = (0), alH > 1 
=>Y.H\Oz{eH - ^ H%Oz{iH)) surjective 

a 

hence by (9) 

H%Opn{i)) H%OziiH)) surjective 

and 

H'+^{OziiH + 1)52)) = (0), aU ^ > 

w{Oz{m) = {Q). 

Now I{Z) is m-regular if and only if W{Xz{m — i)) = (0), i > 1, hence if 
and only if 

H^{Opn{m - 1)) H'^{Oz{m - 1)) surjective 
W{Oz{m-t-l)) = (0), i > 1. 
By the previous remark, this follows provided that 

W+^{Oz{{m -i-l)H- {j + 1)^2)) = (0), if i,j>0,i + j> 1. 
But let us rewrite: 

{m-i-l)H- {j + 1)B2 = Kz+jBi + {m - i - {j + l){d - 1) + r)H 

using (5) and (8). Note that jBi + £H is an ample divisor if £ > l,j > 
0, because \Bi\ is base point free. Therefore by the Kodaira Vanishing 
Theorem, 

H\Oz{Kz+3Bi + m)) = (0), 1,3 > 1, J > 
and provided m = (r + l)((i — 2) + 2, this gives the required vanishing. 
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4 Applications 



From some points of view, the first main problem of algebraic geometry is 
to reduce the study of a general ideal I to that of prime ideals, or the study 
of arbitrary schemes to that of varieties. One way of doing this is to find 
a decomposition of the ideal into primary ideals: i.e. write it as an inter- 
section of primary ideals. But even when non-redundancy is added, this is 
not unique, and usually one actually wants something less: to find its radical 
and perhaps write the radical as an intersection of prime ideals, or to find its 
top dimensional part, or to find its associated prime ideals and their multi- 
plicities. There are really four computational problems involved here which 
should be treated separately: (i) eliminating the multiplicities in the ideal /, 
(ii) separating the pieces of different dimension, (iii) "factoring" the pieces 
of each dimension into irreducible components, and finally (iv) describing 
the original multiplicities, either numerically or by a primary ideal. Three 
of these four problems are the direct generalizations of the basic problems 
for factoring a single polynomial: we can eliminate multiple factors, getting 
a square-free polynomial, we can factor this into irreducible pieces and we 
can ask for the multiplicities with which each factor appeared in the original 
polynomial. There is a fifth question which arises when we work, as we al- 
ways must do on a computer, over a non-algebarically closed field k: we can 
ask (v) for an extension field k' of k over which the irreducible components 
break up into absolutely irreducible components. 

Classical algorithms for all of these of these rely heavily on making ex- 
plicit projections of V{I) to lower dimensional projective spaces. This can 
be done either by multi-variable resultants if you want only the set-theoretic 
projection, or by Grobner bases with respect to the lexicographic order or 
an elimination order, to get the full ideal / fl fc[Xo, ■ ■ ■ Recent treat- 

ments of multi- variable resultants can be found in ||(Jan81j ], [Uha91], and a 
recent treatment of the basis method can be found in ||GTZ88|| . There is 
no evidence that either of these is an efficient method, however, and tak- 
ing Grobner bases in the lexicographical order or an elimination order is 
often quite slow, certainly slow in the worst case. The general experience is 
that taking projections can be very time consuming. One reason is that the 
degree of the generators may go up substantially and that sparse defining 
polynomials may be replaced by more or less generic polynomials. A specific 
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example is given by principally polarized abelian varieties of dimension r: 
they are defined by quadratic polynomials in (4^ — l)-space, but their degree 
here (hence the degree of their generic projection to P^+i) is 4^r! | [[Vlum70a| . 
In fact, any variety is defined purely by quadratic relations in a suitable 
embedding [[Mum70b| . 

Instead of using real computational experience, the fundamental method 
in theoretical computer science for analyzing complexity of algorithms is to 
count operations. For algebraic algorithms, the natural measure of complex- 
ity is not the number of bit operations, but the number of field operations, 
addition, subtraction, multiplication and (possibly) division that are used. 
In this sense, any methods that involve taking Grobner bases for any or- 
der on monomials will have a worst-case behavior whose complexity goes up 
with the regularity of the ideal hence will take "double exponential time". 
However, it appears that this worst-case behavior may in fact only concern 
problem (iv) - finding the primary ideals - and that problems (i), (ii) and (iii) 
may be solvable in "single exponential time" . The idea that such algorithms 
should exist for finding V{I) set-theoretically was proposed in the 1984 lec- 
ture on which this article is based, but turned out, in fact, to have been 
already proven by Chistov and Grigoriev, cf. their unpublished 1983 note 
[ |CG83 |. Their line of research led, in some sense, to the work of Brownawell 



and Kollar, showing the single exponential bound (v^) C / for m = d" 
where d = max( degrees of generators of /). 

Based on this work, Giusti and Heintz ||GH91|| give a singly exponentially 
bounded algorithm for computing ideals such that ^(q^) are the irreducible 
components of V{I) (over the ground field k). The method depends on 
computing what is essentially the Chow form of each component, and leads 
to an ideal defining this variety but not its full ideal. In fact, their qj may 
be guaranteed to be prime except for possible embedded components. 

A direct approach to constructing both a/J and the intersection of the top- 
dimensional primary components of /, denoted Top(/), is given in a recent 
paper by Eisenbud, Huneke and Vasconcelos [[EHV92|| . Their construction of 
the radical uses the Jacobian ideals, i.e. the ideals of minors of various sizes 
of the Jacobian matrix of generators of /. This is certainly the most direct 
approach, but, again they have trouble with possible embedded components, 
and must resort to ideal quotients, hence they need a Grobner basis of / 
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in the reverse lexicographic order. They compute Top(/) as the annihilator 
of Ext™'^''^(^)(fc[Xo, ■ ■ ■ , X„]//, A;[Xo, ■ ■ ■ , X„]), which is readily found from a 
full resolution using Grobner bases. Their algorithm appears to be practical 
in some cases of interest, but still has double exponential time worst-case 
behavior. 

It may turn out to be most effective in practice to combine these ideas. 
Often an ideal under study has regularity far smaller than the geometric 
degree of its top dimensional components; projecting these components to 
a hypersurface requires computing in degrees up to the geometric degree, 
which is wasteful. On the other hand, methods such as those in ||EHV92 



work better in low codimensions, if only because there are fewer minors to 
consider in the Jacobian matrix. Thus, projecting an arbitary scheme down 
to low codimension and then switching to direct methods may work best of 
all. 

This still does not settle the issue of the complexity of calculating V^, 
or, for that matter, calculating the full prime ideal of any subvariety of 
codimension greater than one. Chow form type methods give you an effective 
method of defining the set V{I) but only of generating I up to possible 
embedded components. For this reason, the two schools of research, one 
based on the algebra of /, the other based on subsets of have diverged. 
If we knew, as discussed in the previous section, that the regularity of a 
reduced ideal could be bounded singly exponentially, then we could bound 
the degrees of the generators of \/7, and, using Brownawell-Kollar, we could 
determine a/J up to these degrees and get the whole ideal. But without such 
a bound, it is still not clear whether only V{I) and not \/7 can be found in 
worst-case single exponential time. 

Let's look at problem (iii). Assume you have found a reduced equidi- 
mensional /. To study splitting it into irreducible or absolutely irreducible 
pieces, we shall assume initially it is a hypersurface, i.e. / = (/). Compu- 
tationally, there may often be advantages to not projecting a general J to a 
hypersurface, and we will discuss one such approach below. Geometrically, 
there is nothing very natural about irreducible but not absolutely irreducible 
varieties: from the standpoint of their properties, they behave like reducible 
varieties, except that, being conjugate over k, their components have very 
similar properties. If the ground field k gets bigger or smaller, the set of 
absolutely irreducible components gets partitioned in finer or coarser ways 
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into the /c-components. If one has never done any calculations, one would 
therefore be inclined to say - let's extend k as far as needed to split our alge- 
braic set up into absolutely irreducible components. This is a very bad idea! 
Unless this extension k' happens to be something simple like a quadratic or 
cyclotomic extension of k, the splitting field k' is usually gigantic. This is 
what happens if one component of V{I) is defined over an extension field ki 
of k of degree e, and the Galois group of ki/k is the full symmetric group, a 
very common occurence. Then V{I) only splits completely over the Galois 
closure of ki/k and this has degree e!. The moral is: never factor unless you 
have to. 

In fact, unless you need to deal simultaneously with more than one of 
its irreducible components, you can proceed as follows: the function field 
K = k[XQ, ■ ■ ■ , Xn]/ {f) contains as a subfield an isomorphic copy of ki. 
you find that field as an extension ki = k[y]/{p{y)) of k, and solve for the 
equation of one irreducible component fi G ki[XQ, ■ ■ ■ by the formula 
Normfc,/fc(/i) = /. 

Pursuing this point, why should one even factor the defining equation / 
over fc? Factoring, although it takes polynomial time [[LLL82|| , is often very 



slow in real time, and, unless the geometry dictates that the components be 
treated separately, why not leave them alone. In some situations, for instance, 
pD84|| one may have an ideal, module or other algebraic structure defined by 
polynomials or matrices of polynomials over a ground ring D = k[y]/{p{y)), 
where p is a square-free polynomial. Thus D is a direct sum of extension 
fields, but there is no need to factor p or split up D until the calculations 
take different turns with the structures over different pieces of Spec{D). 

The standard methods of factoring in computer algebra all depend on (i) 
writing the polynomial over a ring, finitely generated over Z, and reducing 
modulo a maximal ideal m in that ring, obtaining a polynomial over a finite 
field; and (ii) restricting to a line L, i.e. substituting Xi = aiXo + bi,i > 1 
for all but one variable, obtaining a polynomial in one variable over a finite 
field. This is then factored and then, using Hensel's lemma, one lifts this fac- 
torization modulo higher powers of m and of the linear space L. One then 
checks whether a coarsened version of this factorization works for /. This is 
all really the arithmetic of various small fields. Geometrically, every polyno- 
mial in one variable factors over a suitable extension field and the question 
of counting the absolutely irreducible components of a variety is really more 
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elementary: it is fundamentally topological and not arithmetic. One should, 
therefore, expect there to be direct geometric ways of counting these com- 
ponents and separating them. Assuming J is a reduced, equi-r- dimensional 
ideal, the direct way should be to use Serre duality, computing the cohomol- 
ogy H^{fly(^j^), where C is the subsheaf of the top-dimensional 

dualizing sheaf of V{I) of absolutely regular r-forms. Its dimension will 
be the number of absolutely irreducible components into which V{I) splits. 
Calculating this cohomology involves two things: algebraically resolving the 
ideal I and geometrically resolving the singularities of V{I) far enough to 
work out f^y(7). Classically, when I = (/) was principal, was called its 

ideal of "subadjoint" polynomials. 

There is one case where this is quite elementary and has been carried 
out: this is for plane curves. One can see immediately what is happening by 
remarking that a non-singular plane curve is automatically absolutely irre- 
ducible, hence one should expect that its singularities control its decompo- 
sition into absolutely irreducible pieces. Indeed, if C C k[XQ, Xi, X2]/{f) is 
the conductor ideal, then is given by the homogeneous ideal C, but with 
degree being shifted to be polynomials of degree d — 3, d the degree of /. To 
calculate H^, assume Xq is not zero at any singularity of V{f) and look at the 
finite-dimensional vector space of all functions k[Xi/XQ, X2/XQ]/ {C + (/)) 
modulo the restrictions g/{XQ~^) for all homogeneous polynomials g of de- 
gree d — 3. This will be canonically the space of functions on the set of 
components of V{f) with sum 0. In particular, it is (0) if and only if V{f) 
is absolutely irreducible. This follows from standard exact sequences and 
duality theory. It was known classically as the Cayley-Bacharach theorem, 
for the special case where V{f) was smooth except for a finite number of 
ordinary double points. It states that V{f) is absolutely irreducible if and 
only if for every double point P, there is a curve of degree d — 3 passing 
through all the double points except P. 

This example gives one instance where a deeper computational analysis of 
varieties requires a computation of its resolution of singularities. We believe 
that there will be many instances where practical problems will require such 
an analysis. In many ways, resolution theorems look quite algorithmic, and, 
for instance, Abhyankar and his school have been approaching the problem 
in this way ||Abh82|] , as have Bierstone and Milman [pM91|| . However, the 



only case of resolution of singularities to be fully analyzed in the sense of 
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computational complexity is that of plane curves. This has been done by 
Teitelbaum [[rei89|| , |[rei9CI|| . His analysis is notable in various ways: he is ex- 
tremely careful about not making unnecessary factorizations, let alone taking 
unnecessary field extensions, and uses the "D" formalism discussed above. 
He describes his algorithm so precisely that it would be trivial to convert it 
to code and, as a result, he gives excellent bounds on its complexity. 
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